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Abstract 

We show that ionization and dissociation can be influenced separately in a molecule with ap- 
propriate external noise. Specifically we investigate the hydrogen molecular ion under a stochastic 
force quantum mechanically beyond the Born-Oppenheimer approximation. We find that up to 
30% of dissociation without ionization can be achieved by suitably tuning the forcing parameters. 

PACS numbers: 33. 80. Gj, 34. 10. x, 03. 65. Sq 
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I. INTRODUCTION 



The anharmonicity of molecular vibrations makes the dissociation of a molecule by irra- 
diation of laser light a relatively difficult task Consequently, high intensity is required for 
dissociation, for instance, I > lO^^W/cm^ for HF in monochromatic fields In half cycle 
pulse fields the situation is improving due to the broad frequency content of this radiation. 
Even in that case, however, intensities / > 10^'^W/cm^ are required to achieve on the order 

n 

of 15% dissociation probability [3|. At such high intensities, the ionization process dominates 
and masks vibrational excitation and dissociation. Chelkowsky et al. ^| suggested that the 
dissociation threshold of a diatomic molecule can be lowered by two orders of magnitude 
using a down-chirped laser, tracing the decrease in spacing of the excited energies levels, 
and hence dissociation without ionization should be possible. In a similar spirit, circularly 
chirped pulses have been used by Kim et al. js] for the dissociation of diatomic molecules. 
They found that the threshold laser intensity is substantially reduced, which may allow to 
achieve dissociation without ionization. 

Recently, dissociation of diatomic molecules under a stochastic force was investigated rep- 
resenting the diatom as a Morse oscillator 6|. A phenomenon, akin to stochastic resonance 
has also been identified in the interaction of an atom with a strong laser field and additional 
noise [3J. 

In the present study, the external force is a sequence of pulses (kicks) at random times, 
each kick carrying an independent weight 0, O]. Both, the weights Fj and the intervals 
between kicks Ati, follow a Poisson distribution with averages 7 and (At), respectively. 
This type of force, similar to white shot noise, has also been used to model the passage of 



ions through carbon foils [9j. Through its inherent time-structure characterized by (At), it 
can couple selectively to bound degrees of freedom with comparable time scales of the noise. 
Given the difficult dissociation without ionization as described above, it is natural to ask if 
stochastic forcing may allow one to separate effectively ionization from dissociation. 

In order to answer this question we have to include both, nuclear as well as electronic 
degrees of freedom in our description. The simplest molecule displaying ionization and 
dissociation is the hydrogen molecular ion Hg" in an external field. Using a two-level approx- 



imation for the electron [10| or a soft-core potential for the Coulomb interaction 
this system has been investigated previously under strong laser fields. 
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As in [13| we refrain from these approximations and investigate this molecule under the 
influence of a stochastic field with the only restriction that the random force is directed 
along the molecular axis. This leads to cylindrical symmetry of the problem and therefore 3 
degrees of freedom (2 for the electron and one for the nuclei) have to be taken into account 
explicitly. On the same level, the interaction of with laser light has been investigated 
recently and a sensitive dependence of the dissociation on the carrier envelope phase has 
been found 1^. Yet, the maximal dissociation probability reported is below 20 %. 

We begin the main text of this paper in section 2 with a description of the noise source and 
a brief explanation, how quantum evolution under a stochastic force is solved for in time. In 
section 3 we discuss the influence of the same stochastic force on an atom (electronic degrees 
of freedom) and a molecule (vibrational degrees of freedom) separately, since the effect of 
coupling strong noise to quantum systems is not widely known. Then in section 4 we describe 
how we handle H.2 numerically, in particular how we extract ionization and dissocication 
from the numerically obtained time-dependent wavefunction. In section 5, results for the 
dissociation and the ionization probability of H2 will be presented and discussed, while 
section 6 concludes the paper. 

II. DESCRIPTION OF THE NOISE SOURCE 



The stochastic force F{t) we consider here is given by 

Nt 

Fit)=Y^^,6it-U), (1) 

i=l 

and stands for a series of random impulses of strength 7j at times ti, i. e., F{t) is a kind 
of white shot noise [l^ responsible for multiple 5— kicks undergone by the molecule, where 
Nt is the number of kicks up to time t controlled by the Poisson counting process Nt. It 
is characterized by the average kicking interval (At) = about which the actual interval 
Ati = ti — are exponentially distributed, similarly as the actual kicking strengths 7j 
about their mean 7, 

P(Ati) = Aexp(-AAti), P{li) = 7"' exp(-7i/7) . (2) 



In analogy to periodically applied half cycle pulses [17|, we restrict our analysis to positive 



7i and assume that % and ti are mutually uncorrelated random variables generated by the 



distributions functions of Eq. ([2]). The determination of F{t) reduces to the construction 



of a stochastic sequence {jijii) which can be done assuming that 
a Poisson sequence of points leading to a delta correlated process 
that the stochastic force constructed has the properties 



;he random times ti form 
. It is easy to show 



{Fit)) = 7A 

{F{t)Fis)) = 27^A5(^-.)+7'A^ (3) 

where () is understood as an average over a sufficient large number of deterministic re- 
alizations j = 1,...N of F{t) in terms of specified stochastic sequences {'~fl''\t[''^). The 
corresponding power spectrum, i. e., the Fourier transform of {F{t)F{s)) , is given by 

Siiu) = 4^ + -f^X''V2^6{Lo). (4) 
V 27r 



These properties reveal the difference between the present stochastic force (white shot noise) 
and a pure white noise which is delta-correlated with zero mean. 

Note that (At) and 7 are the two relevant parameters characterizing the present noise 
source. 

The determination of the time evolution of a system with a deterministic Hamiltonian Hq 
and an additional stochastic driving H = Hq + xF(t) is straight forward by solving first the 
conventional deterministic time-dependent Schrodinger equation for each realization F^^^ of 
the stochastic force, forming from it the desired observable O*^'-* and finally averaging over 
the reahzations, O = N'^ J2f=i . 



III. IONIZATION AND DISSOCIATION UNDER STOCHASTIC FORCES 

To demonstrate that coupling white shot noise to a bound quantum system is sensitive 
to the time scale of bound motion, we briefiy describe how an atom and a diatomic molecule 
responds to noise. To this end we show in figure [1] how the electronic ground state of a 
one dimensional soft-core model of the H-atom and the vibrational ground state of a Morse 
oscillator describing a HF molecule are depopulated under stochastic forcing. The potential 
in the atom case is given in atomic units by 

V{x, t) = Vc^ xF^it) = - J. + xFa(t), (5) 

va;^ + a 
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time [fs] 

FIG. 1: The population of the (electronic) ground state of the hydrogen H atom versus the pop- 
ulation of the vibronic ground state of a Morse oscillator, describing the HF molecule as func- 
tion of time for the same amount of noise. For HF model, /xq = 3.54076 a.u., a = 1.1741 a.u., 
D = 0.225 a.u.M. 



with a = 2 such that the ground state energy corresponds to that of the 3d hydrogen atom. 
In the molecular case it is given by 

V{x) = Vm + f^oxF^it) =D{1- exp{-ax})2 + fioxF^{t) (6) 

with the dissociation energy D, the length scale a, and the dipole gradient fiQ. Note that 



the eigenstates and the eigenenergies of this Morse potential model are exactly known |18 |. 
We expose both systems to the same stochastic external perturbation, as described above, 
by adjusting the molecular force such that fiQF^{t) = Fg^{t). The average time between kicks 
is chosen to equal the electronic period in the hydrogen atom, which in atomic units is given 
by Te = 27r. 

Numerically, we find a very promising indication that the electron in the atom case does 
react much less to the force than the binding of the molecule. This is shown in figure [1], 
where the auto-correlation function of the ground state electron wavefunction for the atomic, 
as well as of the vibrational ground state in the molecular case is shown. While the atom 
is still dominantly in its ground state, the ground state of the molecule has already been 
depopulated by more than 50%. 
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FIG. 2: Schematic ionization and dissociation (including fragmentation) probabilities as function 
of the kick spacing (At)/Te at a given finite time and for a given kicking strength 7. 

In general, we expect for a system S with a characteristic period T| of bound motion 
that noise will couple energy into the system if its kick spacing is smaller than the bond 
period, (At) < Tg. In the opposite case, (At) ^ T|, the bound motion will be insensitive to 
the noise. Since the molecular vibrational period is much larger than the electronic period 
in an atom, T^j^ ^ T^f, there is a window in (At), where dissociation without substantial 
ionization should be possible (see the sketch, fig. [2]). 

We double check this conjecture in the following with the hydrogen molecular ion H2 in 
an external field which is the simplest molecule displaying both, ionization and dissociation. 
Since this problem is sensitive to different timescales, we avoid the Born Oppenheimer 
approximation which could be an (artifical) source of physical effects to be predicted. 




IV. MOLECULAR IONIZATION AND DISSOCIATION DYNAMICS IN H+ UN- 
DER EXTERNAL DRIVING 



We solve the time-dependent Schrodinger equation for the hydrogen molecular ion in a 



linearly polarized laser field neglecting the center of mass and the rotational motion 13|. 



Adding a dipole coupled stochastic force F{t) the time-dependent stochastic Schrodinger 



6 



equation is given by (in atomic units) 

p, R, t) = [i^^ib + K + V + zKFit)] ^(z, p, i?, t) (7) 
with the kinetic energy i^vib = —M'^&^/dR^ of the protons, the electronic kinetic energy 

2dz^ '2dp^ 2pdp 

and the potential energy 

V{p, z, R) = -[p^ + {z- R/2f]-''^ - + {z + i?/2)2]-V2 + II R , (9) 

where (3 = 1/2 + l/(4Mp), k = I + l/(2Mp + 1), Mp is the proton mass in units of the 
electron mass, R is the internuclear distance and z and p are cylindrical coordinates of the 
electron. 

Due to the cylindrical symmetry of the stochastically driven system an expansion of the 
wavefunction in a Bessel-Fourier series in p is performed. The singularities of the Coulomb 
potential is retained, without any softening of the potential. The time-dependent stochastic 
Schrodinger equation (JTj) can then be solved by the standard split operator FFT technique 
{2^. We start from the ground state as initial state \l/(p, 2;, i?, t = 0), represented very 
accurately by the product of the vibrational ground state (f)o{R) and the Isag electronic 
wavefunction of 

The observables we are interested in are the ionization probability 

Pi{t) = l- fi{R,t)dR (10) 
Jo 

and the dissociation probability which is defined without ionization according to 



I'm ax 



PB{t) = fi{R,t)dR. (11) 

Here, /i is obtained by integrating |^'(p, z, -R, t) P over electronic coordinates inside a cylinder 
of radius po = 8 a.u. and height 2x zj + R a.s sketched in figure O The integration area would 
split into two cylinders if i? > 2zi. The splitting, however, is never reached in the present 
calculation due to the smallness of our i?-grid. For z > zj + R/2, the system is considered 
to be ionized since the electron is sufficiently distant to each of the nuclei. The nuclear 
separation defining the onset of dissociation into the H+H"*" channel is taken as Rn = 9.5 
a.u. and [l3|. In this configuration, /i is explicitely given by 



z 




Azi+R/2) npQ 

MR,t) = 2n dz dpp\^iR,z,p,t)\^ 

J-{zi+R/2) Jo 

In our calculations, the electronic grid for z contains 1024 points and extends from 
—50 a.u. to 50 a.u.. A quadratic imaginary potential has been used to avoid unphysical 
reflections at the grid boundaries. The onset of ionization was defined to be at zj = 32 a.u.. 
The nuclear grid consists of 256 points, extending from -Rmin = 0.38 a.u. to -Rmax = 24 a.u. 
and 16 basis functions were used in the Bessel-Fourier expansion for p. 

V. NUMERICAL RESULTS 

Dissociation and ionization probabilities of H2 as described in the previous section have 
been obtained for different parameters of the stochastic forcing. Averaging over around 20 
realizations of the noise was enough to converge the results for all the cases presented below. 
The maximum propagation times for the wave function were always well below 500 fs around 
which which rotational motion of the molecule (not included in the present approach) would 
come into play. 

In figure H] we present results for a fixed strength 7 = 0.9 a.u. of the noise. The average 
spacing of the kicks was varied to cover the switching of the djTiamics from pure ionization 
{{At) /Te ^ 1, figure UK) to pure dissociation {{At) /Te ~ 1, figure Hb). While the ionization 
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FIG. 4: Ionization and dissociation probabilities as a function of time for different average kick 
spacings., with (At)/Te as indicated in panels and for a constant kicking strength 7 = 0.9 a.u. 



probability reaches unity the dissociation probability has a maximum of 25 % as can be seen 
in Figure Ht. 

In Figure [5] the spacing between the kicks is fixed at the intermediate value of (At)/Te = 
0.2, where both, dissociation and ionization, coexist for a certain amount of time. From a) 
to c) we reduce the strength of the external forcing. Reducing 7 reduces the ionization rate 
and correspondingly shifts the peak in the dissociation without ionization to later times. 

What are close to optimal parameters for dissociation without ionization for the stochas- 
tically driven hydrogen molecular ion? This question is answered in figure [61 Here, at 
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time [fs] 

FIG. 5: Ionization (black) and dissociation (red) probabilities as a function of time for different 
noise strength and for fixed spacing {At) /Tf. = 0.2: ionization (a) and dissociation (a'): 7 = 0.9 a.u. 
ionization (b) and dissociation (b'): 7 = 0.75 a.u. ionization (c) and dissociation (c'): 7 = 0.5 a.u. 

fixed final time of 100 fs, the probabilities as a function of the average kick intervals are 
displayed. Beyond {At) /T^ ~ 1 the ionization quickly goes to zero. However, also the dis- 
sociation approaches zero, as the system does not respond to rare kicks. An intermediate 
value of {At) /Te ^ 0.5 turns out to maximize dissociation without fragmentation of the 
whole system. The maximal probability is close to 30 % and therefore higher than any 
comparable value we have found in the literature. We did many additional calculations for 
different final times, as well as for different field strengths and found results qualitatively 
similar to those presented in Figure [6l 

These results confirm the expectations from figure [2] with the important difference, that 
for the real molecule, the onset of ionization limits the maximum dissociation which can be 
achieved. The actual value of the maximum, of course, can only be calculated as we have 
done (see figure [6]). 
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FIG. 6: Ionization and dissociation probabilities as a function of kick spacings at a given finite 
time for a kicking strength of 7 = 0.5 a.u.. 

VI. DISCUSSION AND CONCLUSIONS 

Achieving efficient dissociation in small molecular systems without ionization has been 
a longstanding goal. Here, we have explored an alternative to the standard tools, such as 
chirped laser pulses: the application of stochastic driving in the form of white shot noise 
with a characteristic internal timescale. 

For the hydrogen molecular ion we have shown that dissociation and ionization can be 
well separated by a suitable choice of the noise parameters. A major role being played by 
the internal timescale of the noise, the average interval between kicks. At kicking intervals 
equal to the electronic ground state period of the hydrogen atom the driven system switches 
between being ionized to dissociating already for moderate field strengths below one atomic 
unit. 

Stochastic driving with the specified properties can be realized with chaotic light, where 
a short laser pulse with sufficient band width is shaped with random phases and amplitudes 
of each frequency component 
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